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Abstract 

In this paper we formulate a weighted version of minimum problem 
(1.4) on the sphere and we show that, for K < L, if {<Pk}k=i consists 
of the spherical functions with degree less than N we can localize the 
points (£i,... ,£l) on the sphere so that the solution of this problem 
is the simplest possible. This localization is connected to the discrete 
orthogonality of the spherical functions which was proved in (3|. Using 
these points we construct a frame system and a wavelet system on 
the sphere and we study the properties of these systems. For K > L a 
similar construction was made in paper [4j, but in that case the solution 
of the minimum problem (1.4) is not as efficient as it is in our case. 
The analogue of Fejer and de la Valee-Poussin summation methods 
introduced in [3] can be expressed by the frame system introduced in 
this paper. 
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1 Approximation on the sphere 



Let § 2 denote the three dimensional unit sphere, let V be a finite dimensional 



subspace of L 2 (8 2 ) with dimZ/ = K. With respect to a basis {<fk}k=i of V 
any / £ V has a unique representation 

K 



fc=l 



(1.1) 



For 6 e S 2 , 
Z = 1,...,L. 
Denote by 

and 



L let suppose that we know the values of /(&), 

(1.2) 



* := 



£ C 



(1.3) 



\<pi(£l) ••• <pk(£l)J 

Let us consider the approximation problem: 

Find aeC K with ||f - *a|| 2 < ||f - *a|| 2 for all a £ C K 



(1.4) 



or, equivalently, min llf — ^alU, where II o || 2 denotes the normal euclidian 

aGC fc 



norm. For K < L the solution of the approximation problem (jTT4j) can 
be found using the least squares method, and a is solution of the normal 
equations 

(1.5) 



Assuming that the matrix <& has full rank, i.e. rank(<&) 

a= (* H *)~ 1 * H f. 



K we obtain that 
(1.6) 



If number K is large, then the difficulty in computations is the determination 
of (^ H ^)~ 1 , which needs a great number of operations. 

Problem. The question is: for a given system {</?fc} fe=1 how to choose 
so that the computation of (& H &)~ 1 to be not difficult. The 
most simplest case is if for a good choice of (£i, . . . , £l), the matrix ($^$) _1 
is equal by the identity matrix. 

In what follows we will show that, for K < L, if {<fk}^=i consists of 
the spherical functions with degree less than ./V we can localize the points 
(Cii---)Ci) m this way. This localization is connected to the discrete or- 
thogonality of the spherical functions. Using these points we will construct 
a frame system and a wavelet system on the sphere and we will study the 
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properties of these systems. For K > L a similar construction was made 
in paper |4j, but in that case the solution of the minimum problem (1.4) is 
not as efficient as it will be in our case. The analogue of Fejer and de la 
Valee-Poussin summation methods introduced in [3] can be expressed by the 
frame system introduced in this paper. 



2 Spherical harmonics of degree n 

Spherical harmonics play an important role in Fourier analysis. They are 
restrictions to the sphere § 2 of homogeneous polynomials that are solutions 
of the 3 dimensional Laplace equation. 

Let denote by £ = (sin 9 cos tp, sin 9 sin p, cos 9) G § 2 and / : S 2 — > C, /(£) = 
/(sin 9 cos tp, sin 9 sin ip, cos 9) . 

The set of n degrees spherical harmonics is denoted by H n , with 

dim(H n ) = 2n + 1. (2.1) 
For arbitrary n G N the functions 

Y nk {9,cp) := V2n + 1-Plf l (cos0)e ifc¥ ', fc = -n,...,n (2.2) 

Ifcl 

are called the spherical polynomials (functions) of order n, where P 

n are 

the associate Legendre polynomials. The spherical polynomials of order n 
constitute an orthonormal basis for 7i n , where the ortogonality is with re- 
spect to the scalar product induced by the following (continuous) measure 
on the unit sphere 

2-7T 7T 

/ /(e) d£ := ^ y Jf(e,<p) sin 9 d9 dip, (2.3) 
§ 2 oo 

namely 

2tt TV 

{Y n k,Y m i) = J J Y nk (9,ip)Y m i(9, p) sin 9 d9 dp = 5 mn 5 k i 



(where 5 mn is the Kronecker symbol). 

The addition theorem for spherical harmonics is 

n 

Y nk(0Y nk (v) = (2n + 1) P„(£ • 77), e, rj G S 2 , (2.4) 

fc=— n 

where P n is the Legendre polynomial of n-th degree. 
The function K n : § 2 x S 2 -> R, 

K„(e,7?) = (2n + 1) P„(e • t?) (2.5) 
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is called reproducing kernel because it has the following reproducing prop- 
erty: for arbitrary / E Ti n , 

(f,K n (; V )) = f( V ), with 77 GS 2 . (2.6) 

3 Discrete orthogonality of spherical functions 

In paper [3] we proved that it can be constructed a set of points in [0, n] x 
[0, 2tt] and a discrete measure, so that the orthonormality property of spheri- 
cal functions regarding to the scalar product, induced by the discrete measure 
is preserved. In what follows we summarize these results. 
Let denote by Af E (—1,1), k E {1, . . . , N} the roots of Legendre polyno- 
mials Pjv of order N, and for j = 1, . . . , N let 



(x - Af ) . . . (x - Aff.QQc - Af +1 ) ...{x — X%) 
j [X) (Af - Af ) . . . (Af - Af_,)(Af - Af +1 ) . . . (Af - A#) ' 



be the corresponding fundamental polynomials of Legendre interpolation. 
Denote by 



i 



A%:= J Zf(x)dx (l<k<N) (3.2) 
-l 

the corresponding Cristoffel-numbers. Let consider the set of nodal points 

X := L kj = (6 k , <pj) = (arccos Af , : k = TJf, j = O^N 

(3.3) 

and the weights 

A N 

:= pfe)' (3 - 4) 

These weights are positive numbers (see [2]). On the set X we consider the 
following discrete integral 

. N 2N N 2N ^ N 

/ /d MN == x^/m^ko = EE^^'^) 2 (2ivViv (3 - 5) 



X 



k=l j=0 k=l j=0 



Theorem 3.1 ([3]) Lei N E N ; iV > 1, i/ien f/ie /mzie sei o/ normalized 
spherical functions 

{y nfc :S 2 ^C : fcE {-n,...,n}, nE {0, . . . , iV - 1} } 
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form an orthonormal system on the set of modal point X regarding to the 
discrete integral (|3.5I) . i.e., 



YnkY n 'k' d^N = <W<W, n, n < N, k £ {-n, ...,n}, k' E {-n, . . . ,n'} 



(3.6) 



Theorem 3.2 ([3]) For all f E C(§ 2 ), 



x s 2 



4 Approximation on the sphere corresponding to 
points determined by discretisation process 

Let localize £ on the sphere in the following way: 
X' = <£kj = (sin Ok cos ipj, sin 6^ sin ^ , cos#fc) E §> 2 : 



k = l,N, j=0,2iv}. (4.1) 



This set consists of L = N(2N + 1) points. 

Suppose that we can measure the values of / on the set X' , namely we know 

f = Cf(£io),...,/(W)) T 

Let introduce the following notations: 



In 











V 



\ 



E M 



LxL 



\Zhn(£,n2n)/ 



fl = Jivf = (1/ MJv(Clo)/(Clo), \//«aK6v2aO/(6v27v)) j 

Sn*(0 = *WOvWO, ne{0,l,...,iV-l}, fc€{-n,...,n}, 



/ 500(60) ffl-l(£lo) ••• 5iV-liV-l(6o) \ 

500(61) 0i-i(£n) ••• 5iv-iiv-i(6i) 



E M 



LxN 2 



\9oo(£,N2n) 9l - 1 (£iV2]v) ••• 9n - IN - i{£,N2n) / 

Let formulate the following weighted minimum problem: 
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find a so that min ||f*i — &ia\\ = ||fi — «&ia||. The solution a, according 
to (1.6) is 

a = ((I N &) H (In®))- 1 (I n $) H fx. 

We show that in this case ((In®) H (In&)) is the N 2 dimensional identity 
matrix, so the computation of a is the most simplest possible. Indeed, denote 
(I N &) H (I N 4>) = 



( v4(0,0; 0,0) 



A(o,0; N-l,N-l) \ 



£ M 



N 2 xN 2 



\A(N-1,N-1; 0,0) ... A(JV- 1,AT- 1; AT- 1,JV- l) J 



where 



A(n, t, n,k) = >W(£x,o)^n'fc'(£l,oW(£l,o) + - • •+>W6v,2A0^'fcK^2Ar)A t A r (6v,2Ar)- 

From (3.6) and the discrete orthogonality of spherical functions it follows 
that 



A(n,k; n',k') = / Y nk Y n , kl d^ N = <W<W, 



x 



n, n' G {1, . . . , iV - 1} k G {— n, . . . ,n}, fe' G {— n', . ■ 
from this it follows that 



.,n'}, 



(J JV *) H (7 Ar *) 



/l ... 0\ 
1 ... 



G M 



N 2 xN 2 



\0 ... 1/ 
From (1.7) and (3.5) it follows that 



a = (/iv*)^fi 



Yx_x/d/x N 



/ yiv-ijv-i/d^N 



/ 

which means that that the components of a are exactly the discrete Laplace- 
Fourier coefficients of /. In this case the best approximant can be expressed 
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using the reproducing kernels K n (£,r)), namely 



N-1 



f = *ia 



/ E ( K n(£l,0,V)) Vw((,l,0)f{v) d/JN(l) 
X n=0 
N-1 

X n=0 



JV-1 



/ E ( K n(S,N,2N,v)) VVN((,N,2N)f(ri) d/XN^) 
X n=0 



This gives the idea to study the properties of the set of functions 
where 



^•(•,0= E E ^O y ^(-)= E 



i,m ) 



(4.2) 



n=0 fc= 



n=0 



{rrij}^ =1 is a strictly monotone increasing sequence of positive integers, jo so 
that rrij < N, if j < jo. We will show that with this systems we can generate 
a multiresolution decomposition in the space of spherical polynomials of de- 
gree N, they constitute a frame system at every level of the multiresolution, 
and they generate a wavelet decomposition on the sphere. 

5 Multiresolution decomposition and wavelet spaces 
on the sphere using the nodal points determined 

by X 

Definition 5.1 A sequence {l^j}jeN of finite dimensional subspaces o/L 2 (§ 2 ) 
will be called a multiresolution of L 2 (S 2 ) if the following conditions are sat- 
isfied: 

Vj C V j+1 for all j G N, (Ml) 



closure Vj,\ 



L 2 (S 2 ) 



(M2) 



Usually, a definition of multiresolution includes a condition on the intersec- 
tion of the spaces V j . From (IMlj) it follows immediately that 



V 



o- 
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We define the wavelet spaces VU j as the orthogonal complements of V ' j in 
Vj+l, i.e. 

VUj := V j+1 e Vj, (5.1) 

which means that 

//; . i V y VUj. (5.2) 



Lemma 5.1 Let be a multiresolution analysis of L 2 (S 2 ) and for 

j £ N let IAJ 'j be the corresponding wavelet spaces defined by (|5.1I) . With 
IASq := V\ we have 

oo 

I?(S 2 ) = QlA/ j . (5.3) 

j=0 



Defining the operators Rj and Qj to be orthogonal projections Rj : L 2 (S 2 ^ 
'j and Qj : L 2 Wj 



Vj and Qj : L 2 MJ h j G N we have 



R j+1 = R 1 f + J2Q k f (5.4) 



fc=i 

for / G L 2 (§ 2 ) . We will study an approximation process on the sphere based 
on harmonic polynomials. 

Theorem 5.1 Every harmonic polynomial of degree N over the sphere 8 2 
can be written as a sum of spherical harmonics of degree N, i. e. 

N 

n N (§ 2 ) = @H n . (5.5) 

n=0 

Particularly, any function / G L 2 (S> 2 ) can be approximated by spherical har- 
monics in L 2 (S 2 )-sense up to arbitrary precision. From L 2 (S 2 )-orthonormality 
of the space 7i n for n G N we obtain the dimension 

N-l N-l 
dimn"Ar_i(S 2 ) = dimHn = J^(2n + 1) = iV 2 . (5.6) 

n=0 n=0 

In the finite dimensional Hilbert space TL n with the inner product L 2 (§ 2 ) 
any function Y n G TL n can be represented with respect to an orthonormal 
basis {Y nk } k= _ n as a Laplace-Fourier sum 



Y n = (Y n ,Y nk )Y nk . (5.7) 



k=—n 
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Finally, any Y G ILv-i can be written as 

N-l n 

Y (V)=J2 E ( Y ' Y nk)Y nk ( V ). (5i 



71=0 fc= 



Because of the discrete orthonormality of spherical harmonics the Laplace- 
Fourier coefficients a nk ■= (Y,Y nk ) can be computed using the discrete inte- 
gral, namely 



a nk = {Y,Y nk ) x = jY nk Y^ 



x 

N 2N 

E E Y ( e k>, <Pj) ■ Y nk (9 k/ ,ipj)fi N (z k/j ). (5.9) 

k'=l j=0 

This means that if we know the values of Y G Il y-i on the set X then 
we can compute the exact values of continuous Laplace-Fourier coefficients 
of Y, consequently we know the values of Y on the whole sphere. This 
makes possible to develop a construction of weighted scaling functions in 
case K < L, which is analogous to that of J. Prestin and M. Conrad [4] for 
K > L. 

Substituting a Uik in (15. 8p we obtain that 

JV-l n 

Y=J2J2< Y ' Y ^)x Y nk, YeU N ^(S 2 ). (5.10) 

n=0 k =—n 

From (|5.8j) . (I5.10P and using notation (|2.4|l . (|2.6j) we obtain the following 
reproducing formulas in the space IIw-i(!§ 2 )- Any Y G ILjv-i(§ 2 ) can be 
written in following two ways: 

Y( V ) = /Y(.),f^K n (., V )\ = /Y(.),^K n (-, V )\ . (5.11) 

\ n=0 / \ n=0 / x 

Let {nij} ( jL 1 be a strictly monotone increasing sequence of positive integers. 

TOj— 1 

Then the spaces ISj = n mj ._i 

(§ 2 ) = H„, (j G N) satisfy the following 

n=0 

conditions 

1. Vj C Vj+i, J G N 

2. closure Aj£Li || o ||) = L 2 (§ 2 ). 

Thus {^jjj^i is a multiscale decomposition of L 2 (§ 2 ). We will choose from 
V\ the scaling functions 

{-i mi— l;n 
J 7i=0;fc=— 71 
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and the set of nodal points X' on S 2 , corresponding to the set X from section 
2. 

Let consider jo so that rrij < N, if j < jo- For the scale V j we introduce the 
scaling functions 

mj-l n 

&=J2 E y «* ( 5 - 12 ) 

?i=0 k=—n 

and the "weighted" scaling functions from V j 



= E E ^(O^O- ( 5 - 13 ) 
n=0 k=—n 

where ^ G X' are defined by (14~T]) . 

The advantage of this choice is the following: beside the continuous orthonor- 
mality of the spherical functions of degree less then N they have the discrete 
orthonormality property on set X. 

In this way on the j-th scale the number of the weighted scaling functions 
is equal to the number of the points on which we measure the functions, 
namely L = N(2N + 1). From J2~§ and JIl) it follows that 



<Pj(;$n) = E (2n+l)P„(-,G£) = E ( 5 - 14 ) 

n=0 n=0 

are real valued. 

We summarize some properties of the functions ¥?j(">£zm) ' n ^ ne nex ^ result- 



Theorem 5.2 

1. The functions ¥>.?'(■> have the reproducing property 

(f,<Pj(;& = ftilL) for all f G V h j < j 

(f^ j (;& x = mlL)- 

2. It holds = m,, vMrrvZL) = ™l WvAADU = m r 

3. The function <Pj(',£iL) is localized around i.e. 



mm 



: fevj, Z(^) = i}- 



^ We Aave span {o^-, | G X'} = Z/ J; j < j 
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5. The set [yj VN^wJvji; : C/^ G X'} is a tight frame in V 'j, 

j < Jo, i-e. for every function f £ V j we have 



N 2N 



1=1 m=0 



6. 



J ^(£,^)du;(£) = l, 

s 2 

| ^(£,^)du;(£) = l ^<iV. 



x- 

Proof 

1. On the base of l|5.8|) . every f £ Uj = U mj -i(S 2 ) can be written as 

mj— l n 

/fa) = E E </>*w y »*fa). 

n=0 k=—n 

Then 

rry— 1 n 

</,W(',^))= E E </.l r nfc><l r nfc,Vi(-,^)> = 
n=0 k=—n 
m.j—1 n mj—1 n 

E E V> Y ^ E E ( Y nk,Y n > kf )Y n r k r(^J 
n=0 k=—n n'=0 k'=—n 

Using the orthonormality property (|2.3p we obtain that 



m, — 1 



n 



n=0 fc= 



Using the discrete orthonormality (13, 6p in an analogous way we obtain 
that 



ll^(-,6m)l| 2 = <Vi(-,^),ViO,^)> = VjtiLtlL) = 
rrij — l n rrij — l 

E E Yn^DYnki^rn) = E ( 2n + *) ' 6m) 

n=0 k=—n n=0 

rrij—l m,-l 



£(2n + l)P ft (l)= E( 2n + 1 ) 

n=0 n=0 
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3. Let / G Vj with f{tfj = 1. Then 

rrij— 1 n 

1=EE «nfc(/)^nfc(^). (5-15) 
n=0 k=—n 

Applying the Cauchy-Schwarz inequality we obtain that 

m,-l n mi— 1 „, 



!< ( E E Mf)\ 2 )( E E Y nkW Y nk(tfj), (5-16) 
n=0 n=— n n=0 k=—n 

with equality attained for /, for which the vectors 

{«n fc (/)| ' , {Y nk {tL)\ ' 

{ J n=0;k=—n ^ J n=U;fe=— n 

satisfy the follo wing con ditions: there exists a constant a G C such 
that a„fe(/) = aY nk (^ m ) for n = 0, my - 1, k = -n,n. 

Prom 2. we deduce that a = — y, and the corresponding / = -\ipj(-, 

For all / G i^j with = 1, from (15.16P and the Parseval's equation 

it follows that llfll 2 > 



Thus 



mm 



4. Writing / G Z/j as its Fourier-Laplace sum 

rrij— 1 n 

/ = E E a nk{f)Ynk, 
n=0 k=—n 

on the base of (15. 9|) . a n fc can be computed as 

JV 27V 



fc'=l jr=0 



This means that 



"i-i-l n / N 2N 



f(0 = E E EE/(^)^(4)^(4) y ™^) = 

n=0 fc=-n yfc'=lj=0 / 
N 2N 

= </(0. (5-17) 

fc'=l j=0 

Hence is spanned by the functions <Pj(-,£i>j), Ck'j ^ -^'> 3 — Jo < -W- 
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5. For every function / G U j we write 

m,— 1 



j n 

2 



</,/> = E E «nfc(/)</,^>- 



n=0 k=~n 

From (|5T7D it follows that 

m 3 — 1 n AT 2Af 
n=0 fc=— n fc'=l j=0 



Using 1. we obtain that (^'Pj(-,^k'j)' Y nkJ = Y nk (£ klj ) and /(^-) 
/,^(-,^,)),then 



n N 2N 
n=0 fc=— n fc'=l j=0 

EE(^(-.®m^) E E ^(7)^5) = 

fe'=l j=0 n=0 k=—n 
N 2N 

EE</>^(-'^))/(e£)Miv(e£) = 



fc'=l j=0 
A 2AT 

*AT 

AT 27V 



EE V-MMft 

fc'=l j=0 



This means that 



is a tight frame in ZA,- with frame bound A = 1. 
6. From the orthonormality property (|3.6|) and taking into account that 



13 



^oo (0 = 1 we have 

. m,--l „ 



g 2 n=U k=—n 
rrij — 1 n 

E E | dw(ow^K) 



n=0 k=—n g 2 



m, — 1 



E E ^(O^Ofc = I- 
n=0 k=—n 

The proof regarding to the discrete integral J* x t Pj('^iL) 1S analogous. 

We mention that for the special choice mj = j, with the corresponding frame 
system 

the analogue of Fejer and de la Valee-Poussin summation methods introduced 
in [3] can be expressed. 

6 The harmonic polynomial wavelet space 

We define the wavelet space lASj, j 6 N, j < jo — 1 as the direct sum 

mj+i-l 

VUj= H n . (6.1) 

n=mj 

The dimension of VU j is dim W j = dimi^+i — dimiA,- = m| +1 — m|. Let 
consider the wavelets tpj for VU j 

m j+1 -l n 

i>= E E Y ^ (6-2) 

n=m,j k=—n 

and the "weighted" wavelets 

^jUhn) = E E y nfc(^)^nfc G (6-3) 
n=rrij k=—n 

From the addition theorem (12. 4|) it follows that 

mj+i— 1 n 

^•(-,^)= E E( 2n+1 ) p «(-'^)' ( 6 - 4 ) 

n=rn,j k=-n 

consequently ^'(^C^), with £^ S X' is real valued. We summarize the 
properties of ipji'i^hn) in the following theorem. 
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Theorem 6.1 

1. For every f G l^j, j < Jo ~~ 1 6 -X 7 ^ere are t;aM the 
reproducing properties 

2. T/ie following orthogonality properties hold: 

<^-(-,^),^(-,0> = ' (^•(•,^),^(-,O>x = 

pAT pAT p pAT / aJV 

4- 



( ^| = ; = = = min{||/|| : / G H/,-, = 1}, 



V'j (C/ m ! Cj m ) \j m2 j + l ~ 171 

5. 

span{^(-,^), Cim e X'} = M/„ j < j'o - 1, 

/(0 = </,^-(^-))x. /eM/j, £eS 2 , j<jo<N. 

6. T/ie sei 

is a tag/ii frame in W^j, j < Jo ; *- e - / or ewen/ / G W^j, 

AT 2AT 



/=1 m=0 



7. 



/ ^(e-^)d(O = 0, 
/ • d MN (0 = 0, j < Jo < N, rrij > 1. 



Proof 
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1. Every / G VU j can be written as 

f= E E a nk(f)Ynk, a nk = {f,Y r 



n=rrij k=~n 

hence the scalar product it is equal to: 

m j+l-l n 

(f^jUiL))= E E {f,Y nk )(Y nk ^ j {.^ m) ) = 

n=m,j k=-n 
n rn.j+i — 1 n ' 

E E </' y «*> E E (^,w^(^) = 

n=m,j k=—n n'=mj k'=—n' <■ <■ 

m, +1 -1 



j ~t~ l n 

sirrah 



E E </> y n fe >^(^) = /(^ 



n=m.,- fe=-n 

For the discrete scalar product the proof is similar. 

2. From the definition of <pj, tpj and the continuous and discrete orthog- 
onality of spherical functions it follows these orthogonality properties. 
Indeed, for the continuous case we have: 



^■(■,^),^-(-,0> 



E E^^o E E W0iw($)) = 



n=0 k=—n 



E E E E ^(o.yn^co) y^pwO = °- 

n=0 k=—n n'=mj k'=—n' i i n 

The proof for the discrete scalar product is analogous. 



^(,0ll 2 = (^(-,O,^(-,C)> 



y^£U£j = E E y te™) y (^) = 

n=rrij k=~n 
— 1 tHj+i— 1 

E (2n + l)P„(^-^)= E (2n + l) = K +1 ; 



2 2 
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4. Let / G VUj with = 1, then 



m j+1 -l n 

1 = E E <Xnk(f)Y nk (tfj. (6.5) 

n=m,j k=-n 



Applying the Cauchy-Schwarz inequality we obtain: 



i<( E (Ei«^(/)i 2 )( E E y ^(^) y ^te™))' 

n=rrij k=-n n=rrij k=—n 

with equality for / with 

/ 7\ v f cN \ n = mj,...,m j+1 -l 
a nk (f) = aY nk (^J k = -n,...,n 

a G C constant. Then 

m j+ i-l n 



f = a E E 

with a = (m , 1 _ m , ) . If / € VU j with /(^) = 1, then 



2 > II f ||2 



( m j+l ~ m j) 

Thus 

min{||/||: fi M^./i^S 1} 1 - ''' 



5. Writing / G W^j, as its Fourier-Laplace sum 

"b+i-l n 

/ = E E a nkif)Ynk for m j+ i < AT, 

n=nij k=-n 

N 2N 



ank(f) = (f,Y nk ) x = 

fc' = l jr=0 

This means that 

m-j+l-l n / N 2N 



n=mj fc=-n \fc'=l j=0 



AT 2iV /m.j+i-1 n \ . 

EE E E YnkWi^nkio ) f($ j )nN($ j ) = / fm&-)dn 

k'=lj=0 V n=mj fe=-n / ^ 
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from this it follows that 



6. 



span{^(.,f&): ® n zX'} = W i if 



m j+1 -l n 

= (/./)= E E ^(7)(/.^nfc) = 
n=rrij k=-n 

n I N 2N \ 

E E ^ EE/(fe<^>ft y »* = 

n=rrtj fc=- n \fc'=l j=0 / 

m^+i-l ra N 2N 

E E«^)EE foi(-.*fc).rnfc) = 

n=m,j fe=— n fe'=l j=0 " v ' " v ' 

AT 2AT m i+ i-l n 

EE(/.^('^)>^(^) E E <Xnk(f)Yr*(&i) = 
k'=l j=0 n=m,j k=-n 



AT 2AT 



E E </. *>,■(•, $ 

fc'=i i=o 

AT 2N 

EE|(/V^>^^ 



N - 
J- 



fc'=i i=o 
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